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Abstract 
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production at threshold in the standard large- N limit of the 0(iV)-symmetric (j) 4 
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amplitudes while the effect of loops is the renormalization of the coupling constant 
and mass. The form of the solution is due to the fact that the exact amplitude of the 
process 2— m vanishes at n>2 on mass shell when averaged over the 0(iV)-mdices 
of incoming particles for dynamical reasons because of the cancellation between 
diagrams. I discuss some possible applications of large- N amplitudes, in particular, 
for the renormalon problem. 
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1 Introduction 



The problem of calculating amplitudes of multiparticle production at threshold has re- 
cently received a considerable interest [l|]-[[nj. The explicit results for the tree level 
amplitudes in 4 theory demonstrate the factorial grows which is expected due to the 
large amount of identical bosons in the final state. Further investigations of the tree am- 
plitudes H H introduced a nice technique to deal with the problem. A very interesting 
property of the tree amplitudes which has been recently pointed out by Voloshin JlJ is the 
nullification of the on-mass-shell amplitude of the process 2— m for n>4. This nullification 
has been extended fl5f-|I"4j| to more general models and made it possible to calculate the 



amplitude 1— >n at the one-loop level. A dynamical symmetry which may be responsible 



for the nullification has been discussed by Libanov et al. |15 



In the present paper I shall make an attempt to understand which properties of the tree 
and one-loop amplitudes could survive in the full theory which includes all loop diagrams. 
For this purpose I consider the iV-component 4 theory which is exactly solvable in the 
large- iV limit |TB|| when only the bubble diagrams contribute if there is no multiparticle 
production. 

I derive the set of recurrence relations between the amplitudes of multiparticle pro- 
duction at threshold in the large- N limit of the 0(iV)-symmetric (p 4 theory which sums all 
relevant diagrams with an arbitrary number of loops. I reduce these recurrence relations 
to a quantum mechanical problem which turns out to be possible due to the factorization 
at large N. I find an exact solution to the problem using the Gelfand-Dikii representation 
of the diagonal resolvent of the Schrodinger operator. The result is quite similar to the 
tree amplitudes while the effect of loops is the renormalization of the coupling constant 
and mass. The form of the solution is due to the fact that the exact amplitude of the 
process 2^n vanishes for n>2 at large N on mass shell when averaged over the O(N)- 
indices of incoming particles. This nullification occurs for dynamical reasons because of 
the cancellation between diagrams. 

This paper is organized as follows. Sect. 2 is devoted to the definitions and the 
description of the kinematics. In Sect. 3 I derive the set of the recurrence relations and 
transform it to a quantum mechanical problem. In Sect. 4 I present an exact solution 
and demonstrate that it satisfies the set of equations. The exact solution is possible since 
the diagonal resolvent of the Schrodinger operator has a very simple form. I extend the 
method of calculating the diagonal resolvent to the case of a more general potential of the 
Poschl-Teller type. In Sect. 5 I consider some implications of the results and discuss why 
the calculation of large- iV amplitudes for multiparticle production at threshold can be 
interesting for the renormalon problem. Appendix A contains a proof of the uniqueness 
of the solution. 
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Figure 1: The graphic representations of the multiparticle amplitudes a b (n) (a)) and 
D ab (n;p) (b)). 

2 The definitions 



The 0(iV)-symmetric <f) 4 theory is defined in the 4- dimensional Minkowski space by the 
following Lagrangian 



(2.1) 



c = \{d^ b ){d,<p b ) - X(0V 6 ) - ±A(0V^ 2 

where b = 1, . . . ,N and the summation over repeated indices is implied. 

We denote by a\ h (n) the amplitude of production of n on-mass-shell particles with 
the 0(iV)-indices b\ . . .b n at rest by a (virtual) particle with the 0(iV)-index b and the 
energy nm: 

<... 6 » = (6i... &n|0 6 (O)|o). (2.2) 

It is convenient to multiply this amplitude by the vectors £ bi which results in a sym- 
metrization over the 0(iV)-indices of the produced particles and define 



\n) = al... bn {n)^...i K . 



(2.3) 



The amplitude a b (n) is depicted graphically in Fig. |T]a. 

The second quantity which appears in the set of the recurrence relations between the 
amplitudes is the amplitude D% b ,, bn (n;p) of the process when two particles a and b with 
the 4-momenta p + nq and —p, respectively, produce n on-mass-shell particles with the 
0(iV)-indices b\ . . . b n and the equal 4-momenta q = (m, 0) in the rest frame. We define 
again 

D ab (n;p) = DZ b ... bn (n;p)e i ...e n . (2-4) 

The amplitude D ab (n;p) is depicted graphically in Fig. [Tj£>. The definitions (|2.3j ) and 
coincide for N — 1 with those originally introduced by Voloshin [I]. 0] • 
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Figure 2: The recurrence relation for the amplitude a b (n). The analytical formula is 
given by Eq. (fTl]), 



It is easy to estimate at the tree level that 

a 6 (n)~(A£ 2 )^£ 6 (2.5) 

where £ 2 stands for £ a £ a . Since 

A ~ 1 (2.6) 



in the large- N limit fL6 |, we choose 



e ~ N (2.7) 
for all the amplitudes fl2.5| ) to be of the same order in 1/N. 

3 Recurrence relations at large TV 



The condition (|2.6| ) leaves at large- N only the diagrams with the largest possible number 
of sums over internal 0(iV)-indices which propagate along closed loops of diagrams. The 
proper recurrence relations which extend the usual Schwinger-Dyson equations to the 
case when n particles are produced are depicted for a b (n) in Fig. |2|. The continuous lines 
are associated with propagation of the 0(iV)-indices. Each vertex is the sum of three 
possible permutations of the 0(iV)-indices. This is taken into account in Fig. ^| by a 
combinatorial factor. The graphic notations becomes clear if one introduces the auxiliary 
field cr(x) = <ft 2 (x) which propagates in the empty space inside vertices. 

The recurrence relation of Fig. |2| reads analytically as 

„*(„)- Xm-e V S - a a (ni)a b (n 2 )a b (n 3 ) 

ni,n2,n 3 =odd n^.n^.n^. [n-L lj(n 2 lj(n 3 1J 

9 n\ a a (ni) r d A k hh . , , , 

+Xm~ 2 £ ^ m U J J_; / j-^D bb (n 2 ;k). (3.1) 

n —odd, 
no=even 



: ni!n 2 ! (nf -I) J (2vr 




Figure 3: The recurrence relation for the amplitude D aa (n;p). The analytical formula 
is given by Eq. (fH^). 

Notice that only the summed over the 0(iv)-indices quantity D bb (n2', k) enters this recur- 
rence relation. 

It is the property of large TV that the recurrence relation for D aa (n;p) expresses it via 
a b and D bb again. The recurrence relation for D aa (n;p) is depicted in Fig. [3| 

The recurrence relation of Fig. |3| reads analytically as! 

Am' 4 n! a b ( ni )a b (n 2 )D aa (n 3 ;p) 

\ > / I \2 2 2-^i °n.n\- 



(p + nq) 2 - m 2 nu ^ odd) ' (n{ - l){n 2 2 - 1) 

ng —even 

A ^ „ n\ r d A k 



E W^r^T I ——-^D bb {ni\ k)D™(n 2 ;p) . (3.2) 



(p + nq) 2 - m 2 ni ^ even ' 1 2 n x \n 2 \ J (2tt) 



Eqs. ( |3.1| ) and ( |3.2| ) look very similar to the ones [Q, O] for iV = 1 at the one-loop 
level for a(n) when one needs D(n;p) only at the tree level so that the second term on 
the r.h.s. of Eq. ( ^.2| ) can be omitted. 

To rewrite Eqs. ( |3.1| ) and ( |3.2| ) in a more convenient form, let us introduce, following 
the approach of Argyres et al. [|J, the generating functions 

f(r)=<e mr + y-^-e nm V- 2 (3.3) 
^ 3 n!(n 2 -l) 



and 



Sjab i 

D ab ( T ;p) = - -+ y D ab (n;p)—e nmT m n . (3.4) 

l> n=even "•• 



Eqs. ( |3.1| ) and ( |3 . 2| ) can then be rewritten, respectively, as 

f d 2 ^ 



^-m 2 -v(r)U a (r)=0 (3.5) 



2 Here and below the poles should be understood according to the Feynman prescription m 2 — > m 2 — iO. 
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and 

' d 2 o , , 1 „ 1 



^ - uj 2 - «(t) e ^-^(r; p) = (3.6) 



where e is the energy component of p, p = (e,p 



oo = yp 2 + m? (3.7) 

and 

r d 4 h 



v(t) = \^(t) + \J —D b \r;k). (3.8) 

This transformation from Eqs. ( |3TT| ) and ( |3.2| ) to Eqs. ( p^5| ) and ( |3.6j ) is quite similar to 
the one 0, [12| for iV = 1 at the one- loop level. The extension to the O(N) case was 
considered by Brown || at the tree level and by Smith at the one loop level. 

It is easy to understand using Eq. Q3.4p that v(t) which is defined by Eq. ( |3.8| ) is 
nothing but the sum of the matrix elements: 

v(t) = A £(&!.. . 6 n |0 2 (O)|o) e 1 ■ ■ ■ e nmT • (3.9) 

n=0 n - 

This formula as well as Eqs. ( |3.5D , ( |3.6| ) can alternatively be derived by the use of the 
functional technique introduced for this problem by Brown which relates r to the time 
variable t = xq of a (x) by 

r = it. (3.10) 
The fact that the matrix elements rather than operators themselves appear in 



Eqs. ( |3.5| ) and ( p.6| ) is due to the factorization at large iV (see e.g. Ref. HI?]). It is 
the reason why Eqs. ( |3.5| ) and ( p.6| ) form a closed set at large N. 

The next step in the transformation of Eqs. ( |3.5| ), ( |3.6| ) and ( p.8| ) is to pass for D ab 
to the mixed representation — the coordinate in (imaginary) time and momentum for 
space. One defines 

D:\r,r')= [ ^ e ^D ab (r;p). (3.11) 

J ATI 

In order to see that this quantity is indeed associated with the Fourier transform of 

amplitudes in energy, we notice that the insertion of Eq. (|3.4j ) on the r.h.s. of Eq. (|3.11|) 
yields 

r df 00 

D^ b {T,T') = / — V e (^m)r-er' p ) . ( 3 .1 2 ) 

One recognizes now that e + nm is the energy component of the 4-momentum p + nq of 
the incoming particle with the 0(iV)-index a while e is that of p for b. In particular, the 
free propagator in the mixed representation is 

D - (T , T0= ^/*_J_ = ^_L e -^l (forA = ). (3.13) 
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The equations (|3.6|) and ( |3.8|) can be finally rewritten in the mixed representation as 
follows: 

«(r)} ^Dt b (r, r>) = -5(t - r') (3.14) 

and 

v( T ) = A$ 2 (r) + A / ^^-m^r.r). (3.15) 

27T Jm 2 

To understand the meaning of the summed amplitude D^(t,t'), let us note that 
D^(t, t') has in the index space the structure 



Cbcc\ cbcc 

^(r,r') + ^f- 



^ C (r, r') = U 6c - ^ <£(r, r') + (r, r') . (3.16) 



The amplitudes G T and G s are associated, respectively, with the tensor and singlet O(N)- 
states of two incoming particles. The averaged over the 0(iV)-indices of two incoming 
particles amplitude is 



-^(r,r') = ( 1 - -) G u (r, r ) + -^(r,r') (3.17) 



while the generating function for the symmetrized over all the n+2 0(iV)-indices ampli- 
tudes is given by G s \ 

D a J>(T 1 r')^ = GS(T 1 T>). (3.18) 



One sees that the averaged amplitude which enter Eq. (|3.15|) coincides with G T at large 
N. 

Substituting (|3.16|) into Eq. (|3.14 ), we rewrite it at large N as the following equation 
for Gl(r,r'): 

| A _ _ w(T ) J GT( T> /) = _ 5(r _ /) , (3.19) 
while f (r) is related to G^(t, t) by 

V ( T ) = A$ 2 (r) + ^ / Wcu 2 - m 2 G£(r, r) . (3.20) 

27T Jm 2 

Therefore, only the tensor amplitude enters at large N. 

Eqs. ( |3.19|) , (|3.20D together with Eq. ( |3.5|) form the closed set of equations which will 
be solved in the next section. 



4 The exact solution 



To solve the set of equations (|3.5| ), ( |3.19| ) and ( |3.20|) , let us first look at Eq. fl3.19| ) for 



given v(t). This equation determines the Green function of the Schrodinger operator with 
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the potential v{r) while r plays the role of a 1-dimensional coordinate. In other language 



Gf£(r, t') is the matrix element of the resolvent 

GZ(t,t') = (t\ 1 \t') (4.1) 

— D z + ijj 1 + v 

where D stands for d/dr for brevity. One should take then the diagonal matrix element 
of the resolvent, G^(r, r), in order to substitute into Eq. ( |3.20| ) and to determine v(r) 
versus $ 2 . 

The general solution of this problem for arbitrary v is given by the Gelfand-Dikii 



formula 18 



G:(T,T) = RM = t?M (4-2) 

1=0 u 

where the differential polynomials Ri[v] are determined recurrently by 

1 „_i „V 1 



Ri[v] = - l {-D 2 -v-D- 1 vD) ■- (4.3) 
and the inverse operator is 

D~ x v{t) = f dx v(x) . (4.4) 



Eq. ( |4.3|) stems from the fact that R w [v] obeys the third order linear differential equation 

| (|D 3 -Dv- vD) R w [v] = uPDRjp] . (4.5) 

The polynomials Ri[v] depend on v and its derivatives = (D s v). The first few 
polynomials are 

Ro[v] = \, Ri[v] = ~, R 2 [v] = ^(3v 2 - v") , 

R 3 [v] = - — (10v 3 - lOvv" - h{v'f + v {A) ) (4.6) 
64 

while for r-independent v(r) = Vq one has 

Ru[vo] = - , I ■■ (v = const.) (4.7) 

2a/ U! + v 

which agrees with Eq. ( |3.13| ) at r = r'. 



Let us first show how our equations recover the known results [IB| about the large- N 
limit of the 0(A r )-symmetric 4 theory. We put £ a — > in order to suppress the particle 
production. Then the solution for v(r) is r-independent. Using Eq. (|4.7|) and denoting 

Vq = m 2 R — m 2 , (4.8) 

one rewrites Eq. ( p.2(J| ) in the form 

m 2 = m\ - ^ J^ 2 duj ^u 2 - m\ . (4.9) 



This formula exactly coincides with the standard expression for the bare mass, m, via the 
renormalized mass, uir, at large iV [ j"6| . The meaning of this result is very simple: since 
£ a = 0, the only diagrams which are left in the recurrence relations of Fig. are those with 
the bubble insertions to the propagator. They result solely in the mass renormalization. 

It is convenient to perform the mass renormalization in Eqs. ( |3.5| ), ( |3.19| ) and ( |3.20| ) 
eliminating the r-independent part of v. We introduce 



v(t) + m 2 — m 2 R 



after which everything is expressed in terms of tur and vr{t): 



m 



R-v R 



(r)}$ 8 (r) 







and 



where 



vr[t) 



A$ 2 (t) + 



du^oo 2 -m 2 R (RM-± 



CO 



p 2 + m 2 R . 



(4.10) 

(4.11) 
(4.12) 
(4.13) 



We introduce also the renormalized coupling constant, Xr, which is related to the bare 
one, A, by 

11 N r \Ioj 2 - m 



A A, 



ilT 2 Jm 2 



duo 



H 



00* 



in order to rewrite Eq. ( 4.12 ) in the form 

XrN 



v r (t) = \ r $ 2 (t) + ^ / duo uy/u>*-m? R (RM - 

27r Jm l \ 



Ri[vr] 



2oo 



(4.14) 



(4.15) 



The integral over uo on the r.h.s. becomes convergent after the renormalizations. 

Eq. ( [4.14 ) coincides with the standard renormalization of the coupling constant at 
large iV |T9]. The meaning of renormalization is that one chooses the bare quantities, 



m 2 and A, to be dependent on the cut-off according to Eqs. ( |4.8|) and ( 4.14 ) to make the 
renormalized ones, m 2 R and Xr, to be cut-off-independent. 

I found the following exact solution to Eqs. (|4.11|) and (|4.15|) : 

1 

2oo 



A^ 2 (r) 



Vr[t) 



Auo(uo 2 — m 2 R ) 
Ai?$ 2 (r) 



where 



A 



A 



R 



R 



1 + 



8tt 2 



differs from A# by a finite factor and $ a (r) reads 

C a m R e mRT 



1 _ ^r£, p 2m R T 
i g e 

Some comments concerning the exact solution are in order: 



(4.16) 
(4.17) 

(4.18) 
(4.19) 
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i) Eqs. (|4.16 ), ( 4.17 ) recover the large- A?" limit of the results |7j for the tree level 
G^(t,t) and the one-loop level $ a (r). 



ii) The fact that ( 4.16 ) has a pole only at u 2 = m 2 R means the nullification of the tensor 



on- mass-shell amplitudes 2— m for n>2. An analogous property holds [JTOj] at A^ = 1 
for the tree level amplitudes in the case of the sine-Gordon potential and for the 
4 -theory with spontaneously broken symmetry 
iii) The bare coupling A is related to Xr according to Eqs. (|4.14 ) and ([4.18 ) as follows 



(AAO- 1 = (AfiAO" 1 - & (j^ d^" 2 ^ 2 - m| + l) (4-20) 

The additional finite renormalization is familiar from the one-loop result || for 
N = 1. 

iiii) Since Xr becomes infinite at the negative value Xr = —8tt 2 /N, vr(t) vanishes at 
this point when one goes around the pole singularity in ( |4 . 1 9| ) . This is presumably 
associated with the double scaling limit of the 4-dimensional 0(AQ-symmetric 4 
theory |TJ. 

iiiii) The large- A^ amplitude ( f4.19| ) is real which is related to the property of nullification 



discussed in the item ii). This is the difference with the A^ = 1 one- loop result M. 

Let us explicitly demonstrate the property listed in the item ii) which is crucial for 
our solution by the tree diagrams for the process 2— >4. The diagrams which contribute 
to D bb (4:;p) at large A" are depicted in Fig. [| together with the values of the 4-momenta 
(measured in the units of m) of each propagating virtual particle for the given kinemat- 
ics when the incoming particles are on mass shell. Evaluating the contribution of each 
diagram of Fig. [|a, b and c and summing up, one gets 

a) + 6)+ c ) = -i + i + i=0. (4.21) 

I performed also an analogous calculation for the process 2— >6. 

Let us now verify by direct calculations that fl4.16|) , ( [4.17|) and (f4.14|) is the solution 



to Eqs. ( ^4.15| ), ( |4.11| ). The proof of the fact that this solution is unique for the problem 
of multiparticle production at threshold is presented in Appendix A. 

It is instructive first to verify that $ 2 is an eigenvector of the operator on the r.h.s. of 

Eq. 



| (±L> 2 - v R - D^vrD) $ 2 = \D~ l (§L> 3 - Dv R - v R D) $ 2 = D- l Dm 2 R <$> 2 = m 2 R $ 2 . 

(4.22) 

This formula is derived using solely Eq. ( |4. 1 1| ) for an arbitrary vr(t) and the fact that 
$ 2 (t) vanishes for r — ► oo. 



We look now at Eq. Ql.l7| ) as an ansatz with some (yet to be determined) coefficient 



Xr. Then Eq. ( |4.22|) transforms into the following second-order nonlinear differential 
equation for Vr(t): 

v' R -3v 2 R = Am 2 R v R . (4.23) 
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a) b) c) 



Figure 4: The tree diagrams for D bb (A;p) at large N. The components of the 4-momenta 
are given for the propagating virtual particles (in the units of mass) for the 
on-mass-shell kinematics. Only the energy and one space components are 
written since two other components vanish for all particles. The momenta of 
the incoming particles are (2, y/3) and (2, — \/3), respectively while that of the 
produced particles is (1,0). The cancellation of three diagrams is illustrated 
by Eq. (pi]). 



Using Eq. ( |4.6| ) one can rewrite Eq. fl4.23| ) as 

R 2 [v R ]=m 2 R R 1 [v R ]. (4.24) 

This equation is remarkable because we can apply again the same operator as in Eq. (|4.22|) 
to both sides of Eq. ( |4.24|) to get, according to Eq. ( [4.3|) , R3[vr] on the l.h.s and R2[vr] 
on the r.h.s. which is proportional in our case to v R . Applying the operator several times, 
we get 

Ri[vr] = -\m 2 R l - 2 v R (4.25) 



which gives the formula ( f4.16p by using Eq. 

Substituting Eq. ( |4.16| ) into Eq. ( f4.15| ), one finds that it is satisfied for any r providing 
X R and \r are related by Eq. ( |4 . 1 8[ ) . 

It is still left to solve Eq. ((123]) to find an explicit solution for vr(t). In fact Eq. ( |4.23| ) 
is well known and has the solution which coincides (for a < 0) with the profile of the 
one-soliton solution of the Korteweg-de Vries equation at a fixed value of time j^Of : 



/ 2m r \ 8am 2 R e 2mRT 2m 2 R 

v r (t) = 2D 2 log (1 - a e 2mRT ) = ^ ^ = ^ ^ , 4.26 

y ' V > (l-ae 2m «f swh 2 (m R T + \ log a) V ! 

where a is an arbitrary constant which is related to this fixed value of time and Am 2 R plays 
the role of the asymptotic soliton speed. Since vr{t) is related to $ 2 (t) by Eq. ( [4.171) , 
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one gets 



The expression ( |4.27|) is familiar from the tree level results 0, j3|. It satisfies the 
classical equation for the 0(iV)-symmetric 4 theory with the coupling constant Xr: 

{D 2 -m 2 R -X R ^(r)}^(r)=0, (4.28) 

which can be easily obtained substituting ( }4.17|) in Eq. ( f4.11|) . 

The value of a can be fixed comparing with the £^0 limit when the multiparticle 
production is suppressed. In this limit only the diagrams with one insertion of the bubble 
chain contribute to the recurrence relation of Fig. |[ These diagrams lead, as is discussed 
above, to the mass renormalization which is given by Eq. (|4.8|) . In the given case they 
change the coefficient m in front of the exponential in the first term on the r.h.s. of 
Eq. ( |3.3| ) to be m R . Therefore, the expansion of $ a in £ should start from 

$ a (r) = Cm R e mRT + CO{^ 2 ) (4.29) 

which fixes a = Ar£ 2 /8 and results in Eq. ( |4.19| ). 

Our method of finding the diagonal resolvent can be extended to the equation 

[jT2 - - 3**(t)} G w (t, t') = -5(t - r>) , (4.30) 

where the operator differs from the one in Eq. ( |3.19| ) by the extra coefficient 3 in front 
of vr(t). This operator emerges || in the N = 1 case. First we express G^ir, r) via the 
diagonal resolvent of the operator on the l.h.s. of Eq. ( 4.30| ): 

G w (r,r) =R^v r ] (4.31) 

for vr being the solution of Eq. ( [4.23] ) . We then verify that 

\ (|L> 2 - 3v R - ZD^vrD^ v r = m\v R - -v R , 

| [\D 2 - 3v R - SD-'vrD) v 2 r = Am 2 R v 2 R . (4.32) 

The second equation is remarkable because it shows that v 2 R is the eigenvector. Applying 
the operator several times, one proves by induction that 



%N = ~< + (4' - lM<- 2 (4-33) 



which results using Eq. (|4.31|) in the formula 

R ^ 1 = — - 3 ^ (T) 4- (a U) 

" [ Rl 2lu 4u(lu 2 - m 2 R ) ^ 8lu(lu 2 - m 2 R )(u 2 - 4m|) ' 1 } 
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a) b) c) 

Figure 5: Some large- N diagrams for G s (4;p). 



It is easy to extend the results of the previous paragraph to the resolvent Ru>[ s ^ s t v r] 



for s > 2 which also appears in applications p, || [L3[ . One notices that 



(4.35) 

Therefore, the highest power of vr which emerges in R^I ^^ vr] is v r since the second 
term on the r.h.s. vanishes for I = s. 



5 Discussion 

An interesting property of the multiparticle production at threshold in the O(N)- 
symmetric 4 theory at large N is that the amplitude a b (n) forms a closed set of exact 
equations together with the tensor amplitude G T (n;p). The 1— m amplitude a b (n) is real 
in a perfect agreement with unitarity since non-trivial intermediate states are suppressed 
by 1/N (like the wave-function renormalization starts from the order 1/N). This is similar 
to the case of N = 1 with spontaneously broken symmetry where the 1— m amplitude is 
real to all orders ||. As far as G T (n;p) on mass shell is concern, it vanishes for n>2 for 
dynamical reasons and is real for the process 2— >2 due to kinematics. I would say that 
the explicit results for a b {n) and G T (n;p) are not very instructive for these reasons for 
evaluations of the cross section. Analogously, the contribution of classical solutions to the 
partition function is exponentially suppressed by iV at large N. 

The singlet amplitude G s (n;p) reveals, on the contrary, a surprisingly non-trivial 
behavior even at large N. Some of the diagrams which contribute to G s (4;p) are depicted 
in Fig. [5]. The combinatorics is now different from that of the diagrams of Fig. ^ so there is 
no cancellation of tree diagrams for G s (4;p) on mass shell which happens n>6 according 
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a) b) c 



Figure 6: The diagrammatic representation of some diagrams for G s (4;p) at large N 
integrated over d A p. The diagrams b) and c) are of the renormalon type while 

: stands for the bubble chain. 

to the explicit result J7| for G s (4;p) at the tree level. Each of vertices can be "dressed" 
with bubble chains and each of the lines of outgoing particles can be substituted by the 
exact amplitude a b {n) to obtain a diagram with more produced particles. 

One can integrate G s (4;p) over the 4-momentum p to obtain diagrams of the vertex 
type. Some of them are depicted in Fig. ^|. The most interesting are the diagrams of 
Fig. ^6, c which are of the type of renomalons |2l| and should behave to k-th order 



of perturbation theory as k\. It would be very interesting for this reason to calculate 
G s (n;p) exactly at large N. This calculation should be simpler than at N — 1 since the 
factorization is no longer valid at finite N when one should deal with the whole chain of 
the Schwinger-Dyson equations for multipoint Green functions which were analyzed at 



the tree level by Smith |TT 



Another interesting model to investigate is the case of the matrix Higgs field which is 
described by the Lagrangian 

C = \ tr(c^) 2 - ^ W> 2 - ^ W> 3 - £ tr0 4 (5.1) 

where is generically N x N matrix. The model greatly simplifies as iV — > oo at 

fixed AgiV or X^N when only the planar diagrams survive similar to the 't Hooft limit of 
QCD |[22|| . I have checked that the nullification of on-mass-shell amplitudes holds in the 
case of the matrix cubic interaction for 2—>3 at large- N. This looks quite similar to what 
is observed in this paper for the 0(iV)-case. The large- N limit of the matrix 4 theory 
while being simpler than the N = 1 case is, however, quite non-trivial like the large-iV 
limit of QCD and the cross section will be no longer suppressed by 1/N. 
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Appendix A The proof of uniqueness of the solution 



Let us apply to Eq. (|4.15 ) the operator ~ (|-D 2 — v r — D 1 vrD^ which enters the recur- 



rence relation QL3| ). The crucial observation is that $ 2 is an eigenvector of this operator: 

i (1^ _ VR _ D-'vnD) $ 2 = m|$ 2 . (A.l) 



This formula is derived in Eq. (|4.22| ) for an arbitrary vr(t) using only Eq. Q4.11|) and the 
asymptotics of $ 2 (t) for t —>■ oo. Therefore, one obtains from Eq. ( |4.15| ) 

at? r 1 \ 2^2/ n . X R N f , 3 /~2 2 ( v> r i RaM r i[ v r] 1 \ 

-ARa[v R ] = \ R m R $ (r) + ^- 2 - dwwyw* - m\ IR^vr] — 5 — 3 —J. 

(A.2) 

Here we used Eq. (|4.3|) and the fact that R^vr] satisfies the linear equation 



The application of the operator successively I times yields 



-4R l+1 [v R } = A R m|$ 2 (r) + ^ / du u 21 ^ - m\ Lr^[vr} - £ ?M) . (A.3) 

Notice that the integral over u is convergent for any I. Multiplying by i> 2l+2 and summing 
over I we get finally 



-AuR u [vr] + 2 



A R $ 2 fr) 



v 2 — m R 



XrN r ^ \]u 2 -m 2 R / uZR„[v R ]-vZR u [v R ] V 
2tt 2 J mi ^ OJ 2 [ v 2 - u 2 2 



R 



The last expression can be conveniently rewritten as the contour integral 



- AuR v [v R ] + 2 = — + — du ^ ^ — {v2 _ z){z _ ^ (A.5) 
where C\ encircles the poles at z = u 2 , v 2 and infinity leaving outside singularities of 

Now the idea is to solve Eq. ( |A.5| ) for vR v [v R ] as a function of v rather than as a 
functional of vr(t). The analytic properties of ^/zR/^[v R ] as a function of z are known 
for our process of multiparticle production — it may have poles and branch cuts on the real 
axis since the operator in Eq. (|3.19|) is self-adjoint. They may start at z = (n/2) 2 m 2 R which 
is associated with creation of n particles. In addition, y/zR^[v R ] is complex conjugate 
in the upper and lower half-planes since it is real in the interval —m R <z<m R (below any 
thresholds). One can always compress the contour C\ in Eq. (|A.5|) to encircle the poles 
and branch cuts of \/zR^[v R ] which yields 



a p r i i o A^ 2 (r) \ r N f y^ 2 - m R f dx xlm ^R^[v R } 
- AvR v [v R ] + 2 = - jr — ^ + J du ^ J- {v2 _ x){x _ — } (A.6) 
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where the integral is along the support of 1m y/zR^j[v R }. 

Let us now take the imaginary part of Eq. (|A.6|) for real v. One gets 



-UmuR u [v R \ = A H $ (r)7r5(i/ - m fi ) + — - / do; J — -— . (A.7) 

Z7T Co 1 — U) ) 

This linear equation for Im vR v [v R ] has the unique solution 

lmuR u [v R ] = \ R $ 2 n5(v 2 - m 2 R ) (A.8) 

where X R is given by Eq. fl4.18| ). Using the requirement that zaRj,[0] = 1/2, we unambigu- 
ously obtain 

Ru[v R ] = ±- - . X *f {T \, (A.9) 
2v Ah>{v 2 — m R ) 

which coincides with ( |4.1(j| ). 

This completes the proof of the uniqueness of the solution ( (4.16[ ), (|4.17[ ), ( |4.1S| ), Q4.19| ). 
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